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1 Introduction 

The exclusive production of proton-antiproton (pp) pairs in the collision of two quasi-real 
photons can be used to test predictions of QCD. At LEP the photons are emitted by the beam 
electrons^ and the pp pairs are produced in the process e+e~ — > e'''e~77 e+e~pp. 

The application of QCD to exclusive photon-photon reactions is based on the work of 
Brodsky and Lepage [1]. According to their formalism the process is factorized into a non- 
perturbative part, which is the hadronic wave function of the final state, and a perturbative 
part. Calculations based on this ansatz [2,3] use a specific model of the proton's three-quark 
wave function by Chernyak and Zhitnitsky [4] . This calculation yields cross-sections about one 
order of magnitude smaller than the existing experimental results [5-10], for pp centre-of-mass 
energies W greater than 2.5 GeV. 

To model non-pcrturbative effects, the introduction of quark-diquark systems has been 
proposed [11]. Within this model, baryons are viewed as a combination of a quark and a 
diquark rather than a three-quark system. The composite nature of the diquark is taken into 
account by form factors. 

Recent studies [12] have extended the systematic investigation of hard exclusive reactions 
within the quark-diquark model to photon-photon processes [13-16]. In these studies the cross- 
sections have been calculated down to W values of 2.2 GeV below which the quark-diquark 
model is no longer expected to be valid. Most of the experimental data, however, have been 
taken at such low energies. 

The calculations of the integrated cross-section for the process 77 pp in the angular 
range | cos^*| < 0.6 (where 9* is the angle between the proton's momentum and the electron 
beam direction in the pp centre-of-mass system) and for 1^ > 2.5 GeV are in good agreement 
with experimental results [9,10], whereas the pure quark model predicts much smaller cross- 
sections [2, 3]. 

In this paper, we present a measurement of the cross-section for the exclusive process e''"e~ — > 
e+e"pp in the range 2.15 < W < 3.95 GeV, using data taken with the OPAL detector at 
= 183 GeV and 189 GeV at LEP. The integrated luminosities for the two energies are 
62.8 pb"^ and 186.2 pb'^ 

^In this paper positrons are also referred to as electrons. 
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2 The OPAL detector 



The OPAL detector and trigger system are described in detail elsewhere [17]. We briefly 
describe only those features particularly relevant to this analysis. The tracking system for 
charged particles is inside a solenoid that provides a uniform axial magnetic field of 0.435 T. 
The system consists of a silicon micro-vertex detector, a high-resolution vertex drift chamber, 
a large- volume jet chamber and surrounding ^-chambers. The micro-vertex detector surrounds 
the beam pipe covering the angular range of | cos^| < 0.9 and provides tracking information in 
the r-0 and z directions^. The jet chamber records the momentum and energy loss of charged 
particles over 98% of the solid angle. In the range of |cos^^| < 0.73, up to 159 points are 
measured along each track. The energy loss. dE'/dx, of a charged particle in the chamber gas 
is measured from the integrated charges of each hit at both ends of each signal wire with a 
resolution of about 3.5% for isolated tracks with the maximum of 159 points. The 2;-chambers 
are used to improve the track measurement in the z direction. 

The barrel time-of-flight (TOF) scintillation counters are located immediately outside the 
solenoid at a mean radius of 2.36 m, covering the polar angle range | cos 6*1 < 0.82. The outer 
parts of the detector, in the barrel and endcaps, consist of lead-glass electromagnetic calorime- 
try (ECAL) followed by an instrumented iron yoke for hadron calorimetry and four layers of 
external muon chambers. Forward electromagnetic calorimeters complete the acceptance for 
electromagnetically interacting particles down to polar angles of about 24 mrad. 

The trigger signatures required for this analysis are based on a combination of time-of-flight 
and track trigger information. 



3 Kinematics 

The production of proton-antiproton pairs in photon-photon interactions proceeds via the pro- 
cess 



e+(pi) + e-(p2) ^ e+(p;) + e-(p^) + 7(^1) + 7(^2) ^ e+(p;) + e-(p^) + p(?;) + p(g^) (1) 



where denote the four-momenta {% = 1,2). Each of the two incoming electrons emits a 
photon and the flnal state produced by the two coUiding photons consists of one proton (p) 
and one antiproton (p). The four- momentum squared of the two photons is (i = 1, 2) 



Since the electrons are scattered at small angles and they remain undetected, the four-momenta 
squared of each of the two photons are small, i.e. the photons are quasi-real. In this case the 
transverse"^ component of the momentum sum of the proton and the antiproton in the laboratory 
system is expected to be small whereas the longitudinal component of the momentum sum can 

^In the OPAL right-handed coordinate system the z-axis points along the e~ beam direction, and the cc-axis 
points towards the centre of the LEP ring. The polar angle 9 is defined with respect to the z-axis, and the 

azimuthal angle (p with respect to the .x-axis. 

^In this paper transverse momenta are always defined with respect to the z axis. 




(2) 
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be large. The photon-photon centre-of-mass system is generally boosted along the beam axis. 

The larger the boost, the closer the produced (anti) protons are to the beam direction. This 
feature, combined with the typically low mass of the final state pp, and the low efficiency for 
tracking at small angles, leads to significant acceptance losses for these types of events. 



4 Monte Carlo generators 

The e+e^ — > e"''e~pp events are simulated with the PC Monte Carlo generator which has been 
developed to study exclusive photon-photon processes [18]. The PC Monte Carlo generator 
has been expanded for use in this analysis to simulate the kinematics of exclusive baryon- 
antibaryon final states, e"'"e~ — * e"'"e~pp, e"'"e~ e"'"e~AA, and e"'"e~ e"'"e~S~S"'". A total 
of 20 000 e"^e~ — e+e~pp Monte Carlo events have been generated as a control sample for the 
event selection described in the next section. Similarly, for the trigger and detection efficiency 
determination, a total of 360 000 e+e~ e+e'pp Monte Carlo events have been generated 
in the range of 2.15 < W < 3.95 GeV. The background coming from e"'"e~ — > e+e~7r"^7r~ is 
generated with the PC Monte Carlo program, as is the feed-down background from proton- 
antiproton pairs coming from the reaction e"'"e~ e"'"e~AA pp7r+7r^ where the pions are not 
detected (200 000 events). 

The leptonic photon-photon background processes e"'"e~ — > e'^e~ fi~ , e+e" e+e~e"'"e~ 
and c+e~ e'^e~T~^T~ are simulated with the Vermascrcn generator [19]. The KORALZ 
generator [20] is used to simulate the background processes e+e~ — > /^"""/i", and e+e^ — > r+r". 
The e"'"e~ — > e''"e~ background process is simulated with the BHWIDE [21] generator. Table 1 
lists all the generated background Monte Carlo samples. 

Monte Carlo events are generated at y^s^ = 189 GeV only, since the change in acceptance 
between ^/s^ = 183 and 189 GeV is small compared to the statistical uncertainty of the 
measurement. They have been processed through a full simulation of the OPAL detector [22] 
and have been analysed using the same reconstruction algorithms that are used for the data. 



5 Event selection 

The e+e" e+e~pp events are selected by the following set of cuts: 

1. The sum of the energies measured in the barrel and endcap sections of the electromagnetic 
calorimeter must be less than half the beam energy. 

2. Exactly two oppositely charged tracks are required with each track having at least 20 hits 
in the central jet chamber to ensure a reliable determination of the specific energy loss 
dE/dx. The point of closest approach to the interaction point must be less than 1 cm in 
the r4> plane and less than 50 cm in the z direction. 
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3. For each track the polar angle must be in the range |cos^| < 0.75 and the transverse 
momentum p±_ must be larger than 0.4 GeV. These cuts ensure a high trigger efficiency 
and good particle identification. 

4. The invariant mass W of the pp final state must be in the range 2.15 < W < 3.95 GeV. 
The invariant mass is determined from the measured momenta of the two tracks using 
the proton mass. 

5. The events are boosted into the rest system of the measured pp final state. The scattering 
angle of the tracks in this system has to satisfy | cos^*| < 0.6. 

6. All events must fulfil the trigger conditions defined in Section 6. 

7. The large background from other exclusive processes, mainly the production of e"'"e~, 

and tt'^'tt" pairs, is reduced by particle identification using the specific energy loss 
dE/dx in the jet chamber and the energy in the electromagnetic calorimeter. The dE/dx 
probabilities of the tracks must be consistent with the p and p hypothesis. 

- Events where the ratio E/p for each track hes in the range 0.4 < E/p < 1.8 ^ are 
regarded as possible e'*"e~ e+e~e+e~ candidates. These events are rejected if the 
dE/dx probabilities of the two tracks are consistent with the electron hypothesis. 

- Events where the ratio E/p for each track is less than 0.8, as expected for a mini- 
mum ionizing particle, are regarded as possible background from e+e^ e+c^/i+/i~ 
events. This background is reduced by rejecting events where the dE/dx probability 
for both tracks is consistent with the muon hypothesis. This cut is also effective in 
reducing the tt'^tt" background. 

- The dE/dx probability for the proton hypothesis has to be greater than 0.1% for each 
track and it has to be larger than the probabilities for the pion and kaon hypotheses. 
The product of the dE/dx probabilities for both tracks to be (anti) protons has to 
be larger than the product of the dE/dx probabilities for both tracks to be electrons. 

8. Cosmic ray background is eliminated by applying a muon veto [23] . 

9. Exclusive two-particle final states are selected by requiring the transverse component of 
the momentum sum squared of the two tracks, | Z^p±|^, to be smaller than 0.04 GeV^. By 
restricting the maximum value of , this cut also ensures that the interacting photons 
are quasi-real. Therefore no further cut rejecting events with scattered electrons in the 
detector needs to be applied. Fig. 1 shows the | Z^plp distribution obtained after applying 
all cuts except the cut on \J2p±\'^- 

After all cuts 163 data events are selected, 35 events at s/s^ — 183 GeV and 128 events at 
y^s^ = 189 GeV. The distribution of measured dE/dx values versus the particle momentum for 

the selected data events is shown in Fig. 2a. Background from events containing particles other 
than (anti-)protons is negligible due to the good rejection power of the dE/dx cuts. Since no 
event remains after applying the event selection to the background Monte Carlo samples, the 
95% confidence level upper limits for the number of background events expected to contribute 
to the selected data sample are given in Table 1. 

'^E here is the energy of the ECAL cluster associated with the track with momentum p. 
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Since the pp final state is fully reconstructed, the experimental resolution for W (determined 
with Monte Carlo simulation) is better than 1%. The experimental resolution for |cos^^*| is 
about 0.014. Fig. 2b shows the W distribution for data and Monte Carlo signal events after 
the final selection. The Monte Carlo distributions agree well with the data. 



6 Trigger and detection efficiencies 

The e^e^ e+e^pp events contain only two tracks with momenta in the range 0.4 GeV 
to 2 GeV. Special triggers are required to select such low multiplicity events with only low 
momentum particles and the efficiencies of these triggers must be well known. The e"'"e~ — > 
e"'"e~pp events are mainly triggered by a combination of triggers using hits in the time-of-flight 
counters and tracks. The track trigger takes data on the z coordinate of hits from the vertex 
drift chamber, and from three groups of 12 wires in the jet chamber. The selected events must 
satisfy at least one of the following trigger conditions: 

A Two tracks in the barrel region from the track-trigger (TT). This corresponds to an 
angular acceptance of approximately | cos^j < 0.75. 

B A coincidence of at least one barrel track from the track trigger and a. 6 — (p coincidence 
of a track from the track trigger with hits from the time-of-fiight detector (TOF). The 
barrel track and the track forming the 9 — (f) coincidence are not necessarily identical. 
The angular acceptance of the track trigger is approximately | cos^| < 0.75, whereas the 
acceptance of the 9 — (f) coincidence is | cos 6*1 < 0.82. 

Condition A is highly efficient for events within its geometrical acceptance but it triggers on 
two tracks. Condition B is used to trigger on a single track and to measure the trigger efficiency 
in combination with condition A. 

It was checked that the trigger efficiency does not depend on 0. Under this condition 
and assuming the efficiency of each trigger component for each track to be independent, the 
combined event trigger can be written as: 

£trig = £tt + 2£tt-tof£tt(1 — £tt)- (3) 

Here £tt is the efficiency for one track to be triggered by the track trigger and £tt-tof is 
the efficiency for one track to be triggered by the TT-TOF coincidence. The first term in (3) 
gives the efficiency for both tracks to be triggered by the track trigger (condition A) and the 
second term gives the efficiency for the events not triggered by condition A to be triggered by 
condition B for either of the two tracks. 

Data events arc used to calculate the trigger efficiency. The trigger efficiency is determined 
by considering events in which one track detected in one half of the r0 plane (e.g. < < 
180 degrees) satisfies the required track trigger and time-of-flight 9 — (p matching hit while the 
other track in the other half plane (e.g. 180 < (p < 360 degrees) is used to measure the efficiency 
of the track and time-of-flight triggers. 
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The efficiencies measured from the 163 selected events are £tt = (93.7±1.7)% and £tt-tof = 
(58.5 ± 2.3)%, where the uncertainties are statistical only. This yields an overall event trigger 
efficiency of £trig = (94.7 ± 1.5)%, determined from data only. 

To study the efficiency with a large event sample, events from photon-photon processes 
with two tracks in the final state such as e"'"e~ and fi~^fi~ are used together with the pp events 
to determine the track trigger efficiency as a function of p±. This efficiency exceeds 90% for 
p± > 0.4 GeV and the observed dependence of the track trigger efficiency on p± is found to be 
consistent between electrons and muons. 

In a second step, e"'"e~ — * e"'"e~pp Monte Carlo events in the range 2.15 < W < 3.95 GeV 
are used to obtain the trigger efficiencies as a function of W, £trig(W^)- The Monte Carlo 
events are reweighted here according to the trigger efficiency which has been determined as 
a function of the transverse momentum p± using (3). These reweighted Monte Carlo events 
have been used only to determine the trigger efficiency as function of W. Finally the values 
obtained for eTRiG(W^) ai'e normalized to give the overall trigger efficiency ^trig = 94.7% above 
the threshold p± > 0.4 GeV as already determined from data only. The region W < 2.15 GeV 
is excluded from the analysis because the trigger efficiency drops rapidly below 70%. 

The detection efficiency is determined by comparing the number of Monte Carlo events 
passing all cuts with the total number of Monte Carlo events generated within a polar angle 
I cos6'qenI < 0-6 in the pp centre-of-mass system: 

_ da/dW{\cose*\) 

" da/dW^dcOS^^ENl)' 

where | cos 9* \ refers to the reconstructed polar angle in the pp centre-of-mass system. The 
detection efficiency is about 2% at high W and about 4% at low W. To be able to compare the 
measured cross-section with any given model in bins of | cos^^*| and W, the detection efficiency 
has been determined from the signal Monte Carlo in bins of W and | cos 6* \ . 



(4) 



7 Cross-section measurements 



The differential cross-section for the process e+e — > e+e pp is given by 

dV(e+e- e+e-pp) A^ev(W^, I cos^ |) 

dW^d|cOSe*| ~ >Ce+e-£TRIG£DET {W, \ COS^*]) AW^ A| COS^*| ^ ' 

where A^ev is the number of events selected in each (W, \ cos^*|) bin, £trig is the trigger effi- 
ciency, £:det is the detection efficiency, /2e+e- is the measured integrated luminosity, and AW 
and A I cos^*| are the bin widths in W and in | cos^*|. 

The total cross-section (7(77 pp) for a given value of ^/s^ is obtained from the differential 
cross-section d(7(e'^e~ — > e+e~pp)/dl^ using the luminosity function d>C^^/dl^ [24]: 

, , dale'^C e+e~pp) / dC^^ 

^(7T^PP)= ''7 j#- (6) 
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The luminosity function dC^j/dW is calculated by the Galuga program [25]. The resulting 
differential cross-sections for the process 77 — ^ pp in bins of W and | cos 6* \ are then summed 
over I cos^*| to obtain the total cross-section as a function of W for I cos 6** I < 0.6. 



8 Systematic uncertainties 

The following sources of systematic uncertainties have been taken into account (Table 5): 

Luminosity function: The accuracy of the photon-photon luminosity function used in (6) 
has been estimated by taking into account different models such as the p-pole model [26], 
the Equivalent Photon Approximation (EPA) [27], the Generalized Vector Dominance 
Model (GVMD) [25], the luminosity functions given in [28], and in [29,30]. Each of the 
resulting numbers has then been compared with the VDM form-factor model [25] which 
is used for the final result. We take the largest deviation resulting from these comparisons 
as the uncertainty in the luminosity functions which is about ±5%. 

Trigger efficiency: The trigger efficiency model of (3) has been tested by making a compar- 
ison between the measured relative frequencies of the four different sub-combinations, 

(iVxT, A^tt-tof) — (1)1); (2,0), (2,1), (2,2), which can trigger the events, and the pre- 
dicted values, by using fitted efficiencies for TT and TT-TOF (Section 6). Here A^tt 
denotes the number of tracks triggered by the TT and Att-tof the number of tracks 
triggered by the TT-TOF trigger. Table 2 gives the measured and predicted fractions for 
events with W > 2.15 GeV. The ffi of the two efficiencies Ett and £tt-tof yields a 
of 9.6 over 2 degrees of freedom. Although the fit is poor the overall efficiency result is 
consistent with the other determination. An additional systematic uncertainty of 5.0% 
is therefore assigned to the trigger efficiency. This value is obtained by increasing the 
uncertainties to obtain a normalized of one. 

Monte Carlo statistics: The statistical uncertainty on the detection efficiency due to the 
number of simulated Monte Carlo events varies from 4.5% at low W to 6% at high W. 

dE/dx cuts: The systematic uncertainties due to the dE/dx cuts are determined by recalcu- 
lating the Monte Carlo detection efficiency after varying the dE/dx values: 

- The measured values are shifted by ±lcr, where a is the theoretical uncertainty of 
the measurement. 

- The measured values are smeared with a Gaussian distribution of width O.lcr, which 
is the typical systematic uncertainty of each individual dE/dx measurement, as 
found for pions from K° decays [31]. 

The modified dE/dx values are transformed into weights and the new detection efficiency 
is calculated by applying the event selection on the modified Monte Carlo events. The 
systematic uncertainty assigned to each (W, \ cos^*|) bin is the quadratic sum of the full 

deviation of the detection efficiency with smearing and the average absolute deviation 
value of the ±la shifted values with respect to the original values. 

The systematic uncertainties due to the variation of the dE/dx cuts are larger at high 
values of W. They vary between 0.1% for W < 2.6 GeV and up to 5% for W > 2.6 GeV. 
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Residual background: Residual background can come from non-exclusive production of pp 

pairs in processes like e+e" — > e"'"e~pp7r7r if the pions are not detected. The 
cut eliminates most of this background. The data events are almost coplanar, i.e. no 
data event has an acoplanarity of more than 0.262 rad. To estimate the contribution of 
this background, the shape of the acoplanarity distribution for pp pairs in Monte Carlo 
e+e" — > e+e~AA and e+e~ — > e+e^pp events has been fitted to the data. This yields 
an upper limit of 10 events for this background contribution. An additional systematic 
uncertainty of 6% is therefore taken into account. 

Additional uncertainties due to the measured integrated e+e~ luminosity, the track reconstruc- 
tion efficiency and the momentum resolution for the protons are negligible. The total systematic 
uncertainty is obtained by adding all systematic uncertainties in quadrature. 



9 Results and Discussion 



The measured cross-sections in bins of W are given in Table 3. The average (W) in each bin 
has been determined by applying the procedure described in [32] . The measured cross-sections 
(7(77 — > pp) for 2.15 < W < 3.95 GeV and for |cos^*| < 0.6 are compared with the results 
obtained by ARGUS [8], CLEO [9] and VENUS [10] in Fig. 3a and to the resuhs obtained by 
TASSO [5], JADE [6] and TPC/27 [7] in Fig. 3b. The quark-diquark model predictions [12] 
are also shown. Reasonable agreement is found between this measurement and the results 
obtained by other experiments for W > 2.3 GeV. At lower W our measurements agree with the 
measurements by JADE [6] and ARGUS [8], but lie below the results obtained by GLEO [9], 
and VENUS [10]. The cross-section measurements reported here extend towards higher values 
of W than previous results. 

Fig. 3a-b show the measured 77 — > pp cross-section as a function of W together with some 
predictions based on the quark-diquark model [11, 12, 15]. There is good agreement between our 
results and the older quark-diquark model predictions [11, 15]. The most recent calculations [12] 
lie above the data, but within the estimated theoretical uncertainties the predictions are in 
agreement with the measurement. 

An important consequence of the pure quark hard scattering picture is the power law which 
follows from the dimensional counting rules [33, 34]. The dimensional counting rules state that 
an exclusive cross-section at fixed angle has an energy dependence connected with the number 
of hadronic constituents participating in the process under investigation. We expect that for 
asymptotically large W and fixed | cos^*| 

dt ^ ' 

where n = 8 is the number of elementary fields and t = —W^/2{1 — \ cos^*|). The introduction 
of diquarks modifies the power law by decreasing n ton — 6. This power law is compared to the 
data in Fig. 3c with (7(77 — > pp) ~ vr^2{n-3) ^gj^g three values of the exponent n: fixed values 
n — 8, n — 6, and the fitted value n — 7.5 ± 0.8 obtained by taking into account statistical 
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uncertainties only. More data covering a wider range of W would be required to determine the 
exponent n more precisely. 

The measured differential cross-sections d(j(77 pp)/d| cos 6**1 in different W ranges and 
for I cos ^* I < 0.6 are given in Table 4 and in Fig. 4. The differential cross-section in the 
range 2.15 < W < 2.55 GeV lies below the results reported by VENUS [10] and CLEO [9] 
(Fig. 4a). Since the CLEO measurements are given for the lower W range 2.0 < W < 2.5 GeV, 
we rescale their results by a factor 0.635 which is the ratio of the two CLEO total cross-section 
measurements integrated over the W ranges 2.0 < Vl^ < 2.5 GeV and 2.15 < W < 2.55 GeV. 
This leads to a better agreement between the two measurements but the OPAL results are still 
consistently lower. The shapes of the | cos^*| dependence of all measurements are consistent 
apart from the highest | cos^^*| bin. where the OPAL measurement is significantly lower than 
the measurements of the other two experiments. 

In Fig. 4b-c the differential cross-sections d(T(77 — > pp)/d|cos^*| in the W ranges 2.35 < 
W < 2.85 GeV and 2.55 < W < 2.95 GeV are compared to the measurements by TASSO, 
VENUS and CLEO in similar W ranges. The measurements are consistent within the uncer- 
tainties. 

The comparison of the differential cross-section as a function of |cos^*| for 2.55 < W < 
2.95 GeV with the calculation of [12] at = 2.8 GeV for different distribution amplitudes (DA) 
is shown in Fig. 5a. The pure quark model [2,3] and the quark-diquark model predictions lie 
below the data, but the shapes of the curves are consistent with those of the data. 

In Fig. 5b the differential cross-section da(77 pp)/d| cos^^*| is shown versus | cos^^*| for 
2.15 < W < 2.55 GeV. The cross-section decreases at large | cos6'*|; the shape of the angular 
distribution is different from that at higher W values. This indicates that for low W the 
perturbative calculations of [2,3] are not valid. 

Another important consequence of the hard scattering picture is the hadron helicity con- 
servation rule. For each exclusive reaction like 77 — > pp the sum of the two initial helicities 
equals the sum of the two final ones [35]. According to the simplification used in [11], only 
scalar diquarks are considered, and the (anti) proton carries the helicity of the single (anti) 
quark. Neglecting quark masses, quark and antiquark and hence proton and antiproton have 
to be in opposite helicity states. If the (anti) proton is considered as a point-like particle, 
simple QED rules determine the angular dependence of the unpolarized 77 — > pp differential 
cross-section [29]: 

do- (77 pp) (1 + cos^ 9*) 

d|cose*| °^ (1 -cos^e*)' ^ ' 

This expression is compared to the data in two W ranges, 2.55 < W < 2.95 GeV (Fig. 5a) and 
2.15 < W < 2.55 GeV (Fig. 5b). The normalisation in each case is determined by the best fit 
to the data. In the higher W range, the prediction (8) is in agreement with the data within 
the experimental uncertainties. In the lower W range this simple model does not describe the 
data. At low W soft processes such as meson exchange are expected to introduce other partial 
waves, so that the approximations leading to (8) become invalid [36]. 
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10 Conclusions 



The cross-section for the process e"^e~ — > e'^e^pp has been measured in the pp centre-of- 
mass energy range of 2.15 < W < 3.95 GeV using data taken with the OPAL detector at 
^/s^ = 183 and 189 GeV. The measurement extends to shghtly larger values of W than in 
previous measurements. 

The total cross-section (1(77 —>■ pp) as a function of W is obtained from the differential cross- 
section do"(e"'"e~ e~^c~pp)/dW using a luminosity function. For the high pp centre-of-mass 
energies, W > 2.3 GeV, the measured cross-section is in good agreement with other experi- 
mental results [5,7-10]. At lower W the OPAL measurements lie below the results obtained 
by CLEO [9], and VENUS [10], but agree with the JADE [6] and ARGUS [8] measurements. 
The cross-section as a function of W is in agreement with the quark-diquark model predictions 
of [11,12]. 

The power law fit yields an exponent n — 7.5± 0.8 where the uncertainty is statistical only. 
Within this uncertainty, the measurement is not able to distinguish between predictions for the 
proton to interact as a state of three quasi-free quarks or as a quark-diquark system. These 
predictions are based on dimensional counting rules [33, 34]. 

The shape of the differential cross-section da{'y^ — > pp)/d| cos^*| agrees with the results of 
previous experiments in comparable W ranges, apart from in the highest | cos^^*| bin measured 
in the range 2.15 < 14^ < 2.55 GeV. In this low W region contributions from soft processes 
such as meson exchange are expected to complicate the picture by introducing extra partial 
waves, and the shape of the measured differential cross-section d(j(77 pp)/d| cos 6**1 does 
not agree with the simple model that leads to the helicity conservation rule. In the high 
W region, 2.55 < W < 2.95 GeV, the experimental and theoretical differential cross-sections 
d(T(77 pp)/d| cos^*| agree, indicating that the data are consistent with the helicity conser- 
vation rule. 
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Monte Carlo 
background process 


Number of events 
generated 


Luminosity 

(fb-i) 


upper limit 
at 95 % CL 


e+e^ 


c+c-c+c- [19] 


800 000 


1.0 


< 0.75 


e+e^ 


c+c-fx+fx- [19] 


600 000 


1.0 


< 0.75 




e+e-T+T- [19] 


430 000 


1.0 


< 0.76 


e+e" 


e+e-TT+TT- [18] 


149 000 


0.6 


< 1.29 


c+e~ 


r+r- [20] 


80 000 


1.0 


< 0.74 


e+e~ 


^ /i+/i- [20] 


80 000 


10.0 


< 0.07 


e+e~ 


^ e+e- [21] 


600 000 


1.0 


< 0.74 



Table 1: Number of generated events and corresponding integrated luminosities for the different 
background processes. Since no events pass the selection cuts, 95% confidence level (CL) upper 
limits are given for the number of background events expected to contribute to the selected 
data sample. 



trigger class 


data (%) 


model (%) 


(2,1) 


56.5 ±3.9 


45.0 ± 2.5 


(2,0) 


11.2 ±2.5 


16.0 ± 1.7 


(1,1) 


6.8 ±2.0 


7.3 ± 1.7 


(2,2) 


25.5 ±3.4 


31.7 ± 2.5 



Table 2: Measured and predicted relative frequencies of four different sub-combinations, 
(iVTT,-^TOF) = (2,1), (2,0), (1,1) and (2,2), which can trigger the e+e" e+e-pp event for 
events with W > 2.15 GeV. The uncertainties are statistical only. 
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W range 
(UeV) 


(UeV) 


Events 


dfj^e^e — > e^e ppj/aVy 
(pb/Lre V j 


cr(77 ^ pp) 
(nb) 


z.lo-z.zo 


2.20 


52 


5.46 ± 0.76 


2.69 ± 0.39 ± 0.28 


Z.Zb-Z.Sb 


2.30 


33 


2.95 ± 0.51 


1.53 ± 0.27 ± 0.16 






oz 


9 _i_ n 4c; 
z.oo m u.^o 


1 QQ _i_ n 9fi -1- n 1 4 
i.oy Hz u.zu nz u.i^ 


2.45-2.55 


2.50 


20 


1.69 ±0.38 


0.96 ±0.22 ±0.10 


2.55-2.75 


2.65 


18 


1.01 ±0.24 


0.62 ±0.22 ±0.08 


2.75-2.95 


2.85 


6 


0.28 ±0.12 


0.19 ±0.11 ±0.02 


2.95-3.45 


3.14 


1 


0.26 ±0.26 


0.05 ± 0.01 


3.45-3.95 


3.64 


1 


0.28 ±0.28 


0.05 lg:^^±o.oi 



Table 3: Number of events and cross-sections for |cos^*| < 0.6. Statistical and systematic 
uncertainties are also given. A Poisson asymmetric statistical uncertainty, determined within a 
68% confidence interval, has been calculated for the last two bins of W with only one selected 
event. 



cos 9* 


2.15 <W < 2.55 GeV 


2.35 <W < 2.85 GeV 


2.55 <W < 2.95 GeV 




d|cose*| \^^") 




di cos y^'-") 




dcr(77->pp) / 1 \ 
d| cos 0*1 '^""^ 


0.00-0.10 


49 


4.91 ±0.70 ±0.48 


23 


1.87 ±0.39 ±0.18 


6 


0.80 ±0.33 ±0.08 


0.10-0.20 


43 


4.62 ±0.71 ±0.45 


23 


1.86 ±0.39 ±0.19 


4 


0.51 ±0.26 ±0.05 


0.20-0.30 


20 


2.61 ±0.58 ±0.28 


5 


0.50 ±0.23 ±0.05 


2 


0.23 ±0.17 ±0.03 


0.30 - 0.40 


16 


2.24 ±0.56 ±0.23 


8 


0.90 ±0.32 ±0.11 


2 


0.31 ±0.22 ±0.04 


0.40 - 0.50 


6 


1.40 ±0.58 ±0.19 


8 


1.14 ±0.40 ±0.14 


6 


1.03 ±0.42 ±0.13 


0.50 - 0.60 


3 


0.67 ± 0.39 ± 0.09 


6 


1.28 ±0.53 ±0.15 


1 


1.17 ±0.58 ±0.16 



Table 4: Number of data events selected in each bin of |cos^*| and the measured differen- 
tial cross-sections d(j(77 pp)/d| cos6'*| for three W ranges. The statistical and systematic 
uncertainties are also given. 



Source of Systematic 


Systematic 


uncertainties 


uncertainty (%) 


Luminosity Function 


5.0 


Trigger Efficiency 


5.0 


Monte Carlo statistics {W < 2.55 GeV) 


4.5 


{W > 2.55 GeV) 


6.0 


dE/dxcuts (W^< 2.55 GeV) 


0.1 


{W > 2.55 GeV) 


5.0 


Residual Background 


6.0 


Total {W < 2.55 GeV) 


10.3 


Total {W > 2.55 GeV) 


12.1 



Table 5: Estimates of relative systematic uncertainties. The total systematic uncertainty is 
obtained by adding them in quadrature. 
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Figure 1: Distribution of the transverse momentum balance, | Y,P±\'^, for the e^e^ — > e"'"e~pp 
events in data (black points with error bars) compared to the Monte Carlo simulation (his- 
togram). The distributions are obtained after applying all cuts except the | Z^plp cut. The 
arrow indicates the cut value. The error bars arc statistical only. The Monte Carlo distribution 
is normalized to the number of selected data events. 
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Figure 2: a) Distribution of the specific energy loss, dE/dx, versus the particle momentum p for 
the proton and antiproton tracks selected by applying the final selection. The curves indicate 
the expected mean values for different particle species, b) Invariant mass distribution for the 
163 selected e^e^ e^e^pp events. The dots with error bars are the data, the histogram 
denotes Monte Carlo events. Errors bars are statistical only. The Monte Carlo distribution is 
normalized to the number of selected data events. 
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Figure 3: Cross-sections o-(77 pp) as a function of W. The data and the theoretical 
predictions cover a range of | cos^^*| < 0.6. The data points are plotted at the value of (W). 
a,b) The data are compared to other experimental results [7-10] and to the quark-diquark 
model prediction [12]. The error bars include statistical and systematic uncertainties, except 
for TASSO [5] where the uncertainties are statistical only, c) The data are compared to the 
quark-diquark model of [11] (dash-dotted line), and of [12] (solid line), using the standard 
distribution amplitude (DA) with and without neglecting the mass rrip of the proton, and with 
the predictions of the power law with fixed and with fitted exponent n. The inner error bars 
are the statistical uncertainties and the outer error bars are the total uncertainties. 
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Figure 4: Differential cross-sections for 77 — > pp as a function of | cos 9* \ in different ranges of 
W; a, c) compared witli CLEO [9] and VENUS [10] data witli statistical (inner error bars) and 
systematic errors (outer bars), and b) compared with TASSO [5]. The TASSO error bars are 
statistical only. The data points are shghtly displaced for clarity. 
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Figure 5: Measured differential cross-section, dcr(77 — > pp)/d| cos ^*|, with statistical (inner 
bars) and total uncertainties (outer bars) for a) 2.55 < W < 2.95 GeV and b) 2.15 < < 
2.55 GeV. The data are compared with the point-like approximation for the proton (8) scaled 
to fit the data. The other curves show the pure quark model [2], the diquark model of [11] 
with the Dziembowski distribution amplitudes (DZ-DA), and the diquark model of [12] using 
standard and asymptotic distribution amplitudes. 
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